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ABSTRACT 

A  method  is  presented  for  the  analytic  design  of  linear 
control  systems.    Design  is  accomplished  by  the  solution  of 
an  algebraic  set  of  simultaneous  equations  derived  from  the 
design  specifications.    Specifications  considered  include  the 
system  damping  ratio,  undamped  natural  frequency,  bandwidth, 
steady  state  error  coefficients  and  root  sensitivity. 

The  method  is  shownWcfT5e  fa's|,  accurate,  specific  and 
capable  of  working  with  any  number  of  parameters .    In 
addition,  the  analytic  approach  minimizes  dependence  on 
graphical  techniques  and  trial  and  error  solutions. 
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SYMBOLS  AND  ABBREVIATIONS 

CE  Characteristic  equation 

Ess  Steady  state  error 

G(s)  Forward  transfer  function 

kQ  Acceleration  error  constant 

Ka  Accelerometer  gain 

K  Gain 

Kp  or  kp  Position  error  constant 

ICy  or  kv  Velocity  error  constant 

Kt  Tachometer  gain 

Mpw  Reasonance  peak 

Superscript  or  Specified  or  numerical  value 
subscript  "o" 

w  (omega)  Angular  frequency 

wn  Undamped  natural  frequency 

wb  Bandwidth  (  -3  db) 

wr  Rea sonant  frequency 

tf{f  )  Polynomial  in  g 

Q  General  sensitivity  parameter 

s  Complex  variable  of  the  form  cr  +  jw 

S  Sensitivity 

SPq  Sensitivity  of  P  with  respect  to  Q 


Sn  y      component  of  root  sensitivity 

gwn  wn  component  of  root  sensitivity 

S  **T  Total  root  sensitivity 

jr      (zeta)  Damping  ratio 


1 .     Introduction . 

Linear  control  system  design  can  be  defined  as  "a  set  of 
procedures  for  determining  the  value  of  system  parameters 
which  will  cause  a  physical  system  to  operate  in  a  prescribed 
manner.  "    The  set  of  procedures  noted  above  is  as  broad  as 
the  definition.    One  may  choose  to  perform  the  design  effort  in 
the  time  domain  by  applying  modern  state  space  concepts  or  by 
direct  solution  of  the  system  differential  equations.    Alternatively, 
the  complex  frequency  domain  may  be  selected,  employing  any 
of  the  universally  accepted  procedures  such  as  Nichols  and 
Bode  plots,  root  locus,  Mitrovic's  Method,  etc. 

Notwithstanding  the  achievements  which  have  resulted 
from  direct  application  of  the  above  procedures,  it  must  be 
admitted  that  they  have  inherent  limitations.    With  the 
exception  of  simple  systems,  time  domain  analysis  virtually 
demands  the  assistance  of  a  digital  computer.    In  a  sense, 
this  is  a  disadvantage  because  computers  are  neithe  r  available 
to  all  people  at  all  times  nor  is  it  economically  practical  to 
employ  a  computer  for  all  computation.    The  complex  frequency 
schemes  are  considered  imperfect  because  they  are  time 
consuming,  i.e.  ,  essentially  graphical  in  nature  and  often 
based  on  trial  and  error,  and  normally  limited  to  one  or  two 
variables. 


In  view  of  the  above  discussion,  it  appears  that  there  is 
need  for  still  another  design  tool.    However,  in  order  for  a 
new  procedure  to  even  merit  consideration,  it  must  be 
relatively  free  from  the  faults  and  limitations  cited  above. 
Specifically,  the  method  should  be  analytically  simple,  i.e. 
extend  the  range  of  problems  which  can  be  solved  without 
resorting  to  graphical  techniques  or  the  exclusive  use  of  a 
computer.    It  should  also  be  rapid,  capable  of  handling  several 
variables,  specific  and  accurate.    The  remainder  of  this  paper 
is  devoted  to  the  development  and  discussion  of  a  design 
procedure  which  is  believed  to  have  these  attributes. 


2 .      The  Algebraic  Method:    Concept  and  Development 

Concept 

The  proposed  design  procedure  is  based  on  the  premise  that 
system  design  specifications  can  be  transformed  into  algebraic 
functions  of  the  system  parameters.    If  the  specifications  are 
independent,  N  specifications  will  produce  N  equations  which 
simultaneously  constrain  the  system.    Control  system  design 
would  thus  revert  to  the  relatively  simple  problem  of  solving  a 
set  of  simultaneous  algebraic  equations. 

The  problem  of  transforming  specifications  into  equations 
is  not  as  broad  and  difficult  as  one  might  initially  assume. 
The  most  common  and  practical  complex  frequency  specifications 
are:  undamped  natural  frequency  (wn) ,  darning    ratio  (£")/  steady 
state  error  (Ess),  bandwidth  (wb) ,  reasonance  peak  (M      )  and 
reasonant  frequency  (wr) .    In  many  cases,  both  $    and  wn  are 
specified.    This  is  equivalent  to  specifying  the  location  of  a 
complex  pole  pair  and  generally  implies  that  these  poles  should 
dominate  the  system  response.    If  attention  is  restricted  to  this 
situation,  then  the  specifications  of  primary  interest  are         , 
wn,  W]-,  and  Ess  because  M       and  wr  can  be  at  least 
approximately  expressed  in  terms  of  jr  and  w   .    The  remainder 
of  this  section  is  therefore  devoted  to  the  development  of  methods 


which  can  be  used  to  obtain  system  constraint  equations  from 
specifications  f or  y    ,  wn,  w,   and  Egs. 

Root  Location  (£°    and  wn) 

The  coefficient  equations  associated  with  the  specification 
of  a  complex  pole  pair  can  be  developed  in  two  ways .    The  first 
method  is  based  on  the  root -coefficient  relations  of  an 
algebraic  equation;  the  second  is  based  on  Mitrovic's  Method. 

Assume  the  systems'  characteristic  equation  (CE) ,  obtained 
by  setting  the  denominator  of  the  closed  loop  transfer  function 
equal  to  zero,  is  the  polynomial. 


(2-1) 


<2<    S<     -  O 


where  s  is  a  complex  variable  and  defined  as 


—  c*- 


J  (2-2) 

and  the  a^'s  represent  the  real  coefficients  which  are  functions 
of  the  system  parameters. 

Since  equation  (2-1)  is  of  order  n,  there  must  be  n  roots 
which  satisfy  the  equation.    Therefore,  an  equivalent  expression 
for  equation  (2~1)  is 

C~,  (2-3) 

where  R^  represents  the  ith  root  of  the  polynomial.    By  expanding 
equation  (2-3)  and  equating  the  coefficients  of  like  powers  of  s  in 
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equations  (2-1)  and  (2-3) ,  n  equations  can  be  obtained  which 
relate  the  roots  and  coefficient  of  the  original  polynomial. 

If  Jr*  0  anc*  w0  define  the  location  of  the  specified  complex 
pole  pair,  it  is  permissible  to  set 

(2-4) 

Rz  -  F°  *^>  -  i  ^  ft  —  ft*' 

(2-5) 
Substituting  these  values  of  R,  and  R2  into  the  n  root -coefficient 
relations  obtained  previously  gives  n  equations  in  terms  of  the 
remaining  n-2  roots  and  the  coefficients  of  equation  (2-1).    By 
algebraic  manipulation,  it  is  then  possible  to  successively 
eliminate  each  of  the  remaining  n-2  roots  and  finally  obtain 
two  equations  involving  only  the  coefficients.    Since  the 
coefficients  are  explicit  functions  of  the  system  parameters , 
any  combination  of  the  parameters  which  satisfy  these  equations 
will  generate  the  desired  pole  pair. 

Example  2. 1 
Let  the  system  characteristic  equation  be 

S3  -h  qz  s2-    +  G*S    +  ae>     «  0 

Find  the  relations  which  must  exist  between  the  coefficients  in 
order  to  locate  a  pole  pair  withj^  =  1/2  and  wQ  =10. 

Denote  the  three  roots  of  the  CE  as  R]_,  R2  and  R3.    Then 
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and  cr,   m  /?,>?*  +  R^  (&t  *-  Rz ) 

a&   =    Rt  /?*  /?3 
az   =    R,  *- R*.  -hRs 

Setting  f%,    a     *>    +    j  /  0  /  3/y. 

Ra.    =    S~  -  j  to    /  W 

Then  /©£>  (qz-/q)    -  ^  (1) 

/0(Q,-/0)    =Qo  (2) 

/OORs    *■  ^7o  (3\ 

Equations  (1)  and  (2)  are  the  desired  coefficient  relations; 
equation  (3)  provides  information  about  the  third  pole. 

While  the  root-coefficient  method  works  well  with  low 
ordered  systems  (and  even  provides  information  about  the 
remaining  n-2  roots),  it  requires  an  inordinate  amount  of 
computation  when  applied  to  systems  of  higher  order  because 
the  coefficient  relations  must  be  extracted  from  a  set  of  non- 
linear equations.    This  difficulty  can  be  overcome  (at  the  expense 
of  losing  information  about  the  remaining  n-2  roots)  by  using 
Mitrovic's  Method.    With  this  method,  a  kth  ordered  system  can 
be  handled  as  easily  as  one  of  third  order.    The  method  is  fully 
developed  in  (3)  and  (5).     Only  a  brief  review  is  given  here. 

The  argument  begins  with  equations  (2-1)  and  (2-2)  defined 
previously.    Two  independent  equations  in  a^,  wR  and    ;^rca»n  be 
obtained  by  substituting  equation  (2-2)  into  equation  (2-1)  and 
setting  the  real  and  imaginary  portions  of  the  resulting 
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expressions  equal  to  zero.    Through  additional  manipulation  and 
the  introduction  of  a  new  variable,  the  equations  take  the  form 


G<   UJ*  4<  (?)  =P  (2-6) 


*r-o 


AC-O 


ci  K   oj$   cj> .  .     (f)-O 


(2-7) 


where       ^  (?)  =  -  £* F^  (f)  -f  ^  (y)  J  (2-8) 

with    *0  (f)  -  ^     4  tfJ  -  "I      a*"*      4>.K  <W«  -  4<  Cf)   . 

Other  values  of  2^(3^   )  for  various  values  of  y  are  tabulated 

in  Appendix  I . 

Equations  (2-6)  and  (2-7)  are  known  as  the  generalized 
Mitrovic  equations  and  can  be  used  to  obtain  the  coefficient 
relations  required  to  locate  a  pole  pair  as  specified  by  JT 
and  wQ. 

Example  2 . 2 

Work  example  2.1  using  Mitrovic' s  Method. 

For  a  third  order  system,  the  generalized  Mitrovic  equations 
are 

«o4*  (f) + **  oj^4o  cf)  +  **  u>i  4>t  (?)  +  *9  *>i  <j>zcr)~  o 

Substituting  wn  =  wQ  =  10  and  using  values  for  ^(1/2)  from 
Appendix  I,  the  above  equations  simplify  to 

Pi   =     tO     .  7o  =.  /oo  qz  -  iOOO 
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The  results  are  equivalent  to  those  obtained  previously. 

Although  only  a  single  pair  of  complex  poles  is  being 
considered  here,  it  should  be  noted  that  either  method  is 
capable  of  handling  any  number  of  real  or  complex  poles 
(providing,  of  course,  that  the  total  number  considered  is  less 
or  equal  than  the  order  of  the  system) .    For  example,  assume 
two  complex  poles  pairs  and  two  real  roots  of  a  sixth  order 
system  are  specified.    The  root -coefficient  method  will  provide 
six  constraining  equations  by  direct  substitution  of  the  given 
values.    If  Mitrovic's  Method  is  used,  four  constraint  equations 
dependent  oh  the  complex  poles  can  be  obtained  by  repeated 
use  of  equations  (2-6)  and  (2-7).    The  remaining  two  equations 
which  are  dependent  on  the  real  roots  are  produced  by  repeated 
use  of 


C-')*^*<7£    *    O 


<-=o 


(2-8) 


obtained  by  letting  s  =  -gr*  in  equation  (2-1) . 

Steady  State  Error  (Ess) 
There  are  three  different  steady  state  error  constants  which 
may  be  specified.     They  are 

position  error  constant    =  kp  =  lim  G(s)  (2-9) 

s-*-o 
velocity  error  constant    =  k     =  lim  sG(s)  (2-10) 

s-»o 
acceleration  error  constant    =  ka  =  lim  s^G(s)  (2-11) 

s**o 
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The  above  definitions  are  only  valid  for  systems  which  have  all 
their  roots  in  the  left-half  of  the  s-plane  (because  the 
definitions  are  merely  restatement  of  the  final  value  theorem 
for  complex  variables)  and  unity  feedback  around  the  forward 
loop  transfer  function,  G(s).    To  satisfy  the  last  requirement, 
it  may  be  necessary  to  manipulate  the  block  diagram  into  an 
equivalent  G(s)  block  with  the  required  unity  feedback. 

The  constraint  equation  is  obtained  by  equating  the  specified 
value  of  Ess  with  the  algebraic  expression  resulting  from 
application  of  equation  (2-9),  (2-10)  or  (2-11).    The  procedure 
is  demonstrated  in  the  following  examples. 

Example  2.3 

For  the  system  shown  in  figure  2-1,  obtain  an  expression 
between  K,  Z  and  P  f or  k    =  6. 

For  the  given  system 


<U*  (S)    =  «*■*<***> 


52  (*+s)(s+C)(*+R) 

Applying  equation  (2-11) 


Ji  =  .^  **  Q^cs)  = 


.    <z 


which,  for  ka  =  6,  reduces  to:       18 OP  =  KZ 

Example  2.4 
For  the  system  shown  in  figure  2-2,  derive  an  equation  in 
terms  of  K,  ,  K2  and  K3  which  will  give  Kv  =  8 . 
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Since  the  system  does  not  have  unity  feedback,  insert 
fictitious  positive  and  negative  feedback  paths  around  the 
system.     Geq(s)  can  then  be  found  by  including  the  positive 
feedback  branch . 


Sx(H-KtK3)-t-sfa+K,K2.)    -K, 

Applying  equation  (2-10) 


Therefore  the  specification  cannot  be  satisfied. 

Bandwidth  (wij 
Bandwidth  will  be  defined  as  the  angular  frequency,  w,, 
at  which  the  magnitude  of  the  closed  loop  transfer  function  has 
a  value  of  0.707,  i.e.  ,  -  3  db  bandwidth.    The  transformation 
of  the  specification  into  the  desired  constraint  equation  can  be 
accomplished  by  following  the  usual  procedures: 

A.  Find  the  closed  loop  transfer  function  and  set  s  =  jw. 

B.  Find  the  magnitude  of  the  resulting  expression  and 
equate  it  to  0.707. 

C.  Substitute  w  =  w^  and  simplify. 

Example  2.5 
For  the  system  shown  in  figure  2-3,  derive  a  relation  in 
terms  of  K  and  IC  which  assures  the  system  has  a  bandwidth 
of  10  radians/second. 
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Following  the  prescribed  procedure: 


Ccs)      _ 


RCju>) 


iK-co*)*  T  io*(£-+KKt)Z 


For  w  =  wb  =  10 ,  /  O  (juS)  lz    __     II 

\Rtju>)l  /z 


•  • 


<2    =   /0*-ZOOK  +/00(S-+KKt)Z' 


The  constraint  equation  obtained  in  the  preceeding  example 
for  a  second  order  system  is  typical  of  the  expressions  which 
are  developed  from  bandwidth  specifications.    Since  the 
expressions  increase  in  complexity  for  higher  ordered  systems 
or  for  systems  having  closed  loop  zeroes,  most  design  problems 
involving  bandwidth  specifications  will  usually  require  the 
assistance  of  a  computer  or  a  graphical  technique.     (An  excellent 
graphical  procedure  based  on  Mitrovic's  Method  is  given  in  (4).) 
If  such  aids  are  not  available,  a  convenient,  and  relatively  fast, 
alternative  solution  is  to  initially  ignore  the  bandwidth 
specification  and  design  the  system  to  meet  the  remaining 
specifications.    The  results  can  then  be  inserted  into  the  system 
bandwidth  equation  to  determine  if  they  also  satisfy  the  bandwidth 
requirement.    If  the  specification  is  not  satisfied,  the  magnitude 
of  the  deficiency  may  suggest  a  more  reasonable  structure  for  the 
next  attempt. 
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BLOCK  DIAGRAM  FOR  EXAMPLE  2.3 
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BLOCK  DIAGRAM  FOR  EXAMPLE  2  . 4 
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FIGURE  2.3 


BLOCK  DIAGRAM  FOR  EXAMPLE  2  .  5 
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3 .     Application  of  the  Algebraic  Method . 

As  stated  previously,  the  basis  of  the  proposed  design 
method  is  the  transformation  of  N  specifications  into  a  set 
of  N  simultaneous  constraint  equations  that  are  explicit  functions 
of  the  system  parameters.    Since  the  specifications  are 
independent,  a  specific  design  will  generally  require  that  the 
system  have  N  parameters.    The  insertion  of  sufficient 
parameters,  usually  compensating  elements,  is  an  integral  part 
of  the  design  procedure. 

It  should  be  noted  that  the  mere  selection  of  sufficient 
parameters  does  not  guarantee  a  satisfactory  solution.    A 
converse  of  this  statement  is  also  true.    An  unsatisfactory 
solution  does  not  imply  that  a  solution  does  not  exist;  one  may 
possibly  be  obtained  by  employing  a  different  set  of  parameters. 
These  statements  suggest  that  there  is  some  trial  and  error 
associated  with  the  procedure.    This  is  correct.    However,  the 
amount  of  trial  and  error  required  can  be  minimized  through 
experience  and  application  of  control  theory  fundamentals.    In 
this  respect,  root  locus  sketches  are  particularly  valuable  in 
obtaining  clues  to  the  nature  of  the  required  parameters. 

The  algebraic  method  has  the  inherent  capability  to  handle 
any  number  of  specifications  or  parameters.    Therefore  a 
parametric  solution  can  generally  be  obtained  by  simply  inserting 
more  parameters  into  the  system  than  there  are  design 
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specifications.    If  M  parameters  are  selected  for  use  in  the 
design  of  a  system  subject  to  N  specifications,  for  M>N,  the 
results  will  be  parametric  and  expres sable  in  terms  of  any  M  -  N 
parameters. 

Example  3. 1 

Given  G(s)  =  K/s(s  +  4) ,  design*  a  system  which  will 
produce  a  pair  of  complex  poles  defined  by  Jp  Q  =  1/2  and  wQ  =  4 
with  Kv  =  4  . 

It  will  be  stated  without  proof  that  K  =  16  satisfies  all 
specifications. 

This  seemingly  trivial  example  has  major  significance. 
Specifically,  one  parameter  successfully  satisfies  three 
specifications  which  is  a  direct  contradiction  of  the  previous 
discussion.    It  would  appear  that  adherence  to  the  general  rule 
of  N  parameters  for  N  specifications  could  possibly  result  in 
the  incorporation  of  redundant  material  into  the  system.    This 
does  not  seem  likely  because  of  the  following  observation: 

In  all  cases  studied,  parameters  which  did  not  assist  in 

the  solution  tended  to  be  eliminated. 
Examples  3.2  and  3.3  will  be  used  as  illustrations  of  this 
observation. 

Example  3.  2 

Insert  cascade  compensation,  Gc(s)  =  (s+Z)/(s+P),  into 
the  system  given  in  example  3.1  and  redesign  the  system  to 


*Note:    All  designs  will  employ  unity  feedback  between  output  and 
input  (in  addition  to  any  other  specified  feedback) . 
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i e e t    b  !i  e    r,ene    s p© c  1    L c n t  i o  n  °  . 


^69  fs)   *         _K*JL«S 


Applying  equation  (2-10): 


K*  -  -4*  (1) 


^P 


System  CE: 

S3  +&(++/>)    -/-  &  (4P+K  )-*.<*    «  O 

Mitrovic  equations:        ^/Vv-Pj  =   V  /=»  -r  K  (2) 

/CfL    ~  t&(*+P)-m  (3) 

A  simultaneous  solution  of  equations  (1) ,  (2)  and  (3)  gives: 

K=  16  P  =  Z 

Since  the  pdle  and  zero  of  the  filter  are  coincident,  the  filter 
performs  no  function  and  can  therefore  be  eliminated. 

Example  3 . 3 

For  G(s)  =  50/s(s  +  5),  it  is  possible  to  locate  poles  at 
-5  * J  5  (JT  0-  0.707  and  wQ  =  5 (2) l/^2)  by  using  either 
tachometer  feedback  (Kj.  =  0. 1)  or  a  cascade  compensator  (Z  =  5, 
P  =  10) .    In  each  case,  Kv  =  5 .    Design  a  system  using  both 
tachometer  feedback  and  a  series  compensator  to  satisfy  the 
same  specifications. 

For  combined  tachometer  feedback  and  cascade  compensation: 


^(s,    «-  -**<S+*} 


S3    *- 5*  ( P+S  +  SOKr)  +s(5P+SOKt2) 

Applying  equation  (2-10):     <v  =  Soz/(£TP  +  <>OKtI)      (1) 
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System  CE: 

Mitrovic  equations:  So  2  =  SO  (P+£"+  SOKr)"Soo         (2) 

£F>  i-  &OKrlk  +SO  «•  /6  (P-^'S  '+■ £"0  -  T )  "  &C       (3) 

Simplifying:  Z  =  5  (4) 

P  =  10(1  -  5KJ.)  (5) 

The  result  obtained  in  this  example  is  doubly  significant. 
First  of  all,  it  indicates  that  the  procedure  is  capable  of 
detecting  conditional  as  well  as  complete  redundancy.    In 
example  3.2,  the  filter  was  redundant  for  all  values  of  the 
compensator  pole  and  zero.    In  this  example,  however, 
equation  (5)  shows  that  the  filter  is  redundant  only  when 
Kj.  =  0.1.    Similarly,  the  tachometer  may  be  discarded  when  P  =  10. 
The  second  point  of  interest  is  the  fact  that  a  unique  solution 
may  not  exist.    In  this  case,  any  positive  values  of  P  and  Kj. 
which  satisfy  equation  (5)  also  satisfy  the  specifications.    Which 
solution  is  best?    At  this  point,  it  is  a  matter  of  conjecture. 
However,  this  example  will  be  rediscussed  in  Section  4  after  a 
means  for  determining  the  best  design  has  been  developed. 

Example  3.4 

/^ 

Given  the  forward  loop  transfer  function,  G(s)  =— ; — r-z — r 

s  (j^'/i;(s  +  7; 

design  series  compensation  to  give  roots  with  T  0~  1/2  and 
wQ  =  6.    The  minimum  acceptable  value  of  Kv  is  4. 15  . 
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For  cascade  compensation: 


Geo  Cs)  = 


K  (s-hlt) 


S*+  7>Sf>33-*-  S*(7.&P+3.£)  +  3-S~P>s 


Applying  equation  (2-10):      Kv  -  >-f.l£  -      **     —  (1) 

0~  S  f~ 


System  CE: 

Mitrovic  equations:  <it    m  &»+  p-  /  i+o£-  (2) 


K  =    Hh£P—  /9-T         (3) 

Solving  equations  (l) ,  (2)  and  (3)  gives:      K=  1475 

P  =  37.9 
Z  =  0.373 
Example  :3 . 5 
The  bandwidth  for  the  system  in  example  3.4  is  about  6.5 
radians/second.    Redesign  the  system  so  that  the  bandwidth  is 
reduced  to  4  radians/second. 

The  introduction  of  another  specification  requires  the 
insertion  of  another  parameter.    Use  tachometer  feedback  around 
the  original  plant  and  series  compensator.    Then 

Ge  (s)  - gjifrgj 


Applying  equation  (2-10):    /<v  **/./£"—  — — (l) 

3*&P  -t-K.IC.TT. 


System  CE: 

S 4  +  s3  (y.  S-,.P)  t-  Sx  (3 .  S +  ?.Sp-t-  KK-r)  +S (s-SP+KKYt+K)  +-Ki~0 

Mitrovic  equations:     Km  .A*jb.f*><^l*J-^rt    {2) 

m 


System  bandwidth  equation  for  w^  =  4: 


1/2  = ^^- 1—1^2 .       (z/; 

Solving  equations  (1)  through  (4): 

K=  252  P  =  5.93 

Kt  =  0.16  Z  =  1.015 
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4.      Sensitivity  as  a  Specification. 

Attention  is  now  directed  to  the  question  posed  during  the 
discussion  of  example  3.3.    If  more  than  one  design  satisfies 
the  specifications,  which  is  the  best  design?    The  fact  that  a 
specific  solution  was  not  obtained  indicates  the  need  for  another 
equation  which,  for  the  system  being  developed,  implies  the  need 
for  another  specification.    Since  the  design  specifications  have 
already  been  incorporated  into  the  solution,  it  follows  that  the 
additional  specification  should  be  general  in  nature  and 
applicable  to  all  systems. 

One  choice  for  a  general  control  system  specification  is 
sensitivity  which  will  be  defined  as 

C*    =       d  (Jtrn    P)       „  dP/r       .  *     (Sit)  (4-1) 

p 
where  Sq  is  read  as  "the  sensitivity  of  a  function  P  with  respect 

to  a  parameter  Q.  "    (This  definition  is  preferred  by  the  author 

because  it  imparts  more  physical  meaning  than  most  other 

definitions  of  sensitivity  found  in  the  literature,  i.e. ,  it  states 

that  sensitivity  is  essentially  the  ratio  of  the  percentage  change 

in  the  function  to  the  percentage  change  in  the  parameter.)    The 

reason  for  this  choice  is  the  fact  that  a  sensitivity  specification 

will  force  the  final  design  to  be  realistic.    Specifically,  it  requires 

the  consideration  of  changes  in  the  parameters  of  the  physical 

system  which  will  result  because  of  tolerances,  aging  and 

environment.    On  this  basis,  the  best  design  is  the  one  which  not 
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only  satisfies  the  design  specifications,  but  is  also  able  to 
maintain  the  desired  performance  in  spite  of  undesired  parameter 
variations. 

Root  Sensitivity- 
Use  of  equation  (4-1)  requires  definition  of  the  function  P. 
Since  sensitivity  is  being  introduced  as  a  means  of  insuring 
compliance  with  the  design  specifications,  P  should  ideally 
be  a  function  of  all  the  design  specifications.    However,  in 
this  presentation,  only  the  root  specification,  i.e.  ,  j?  Q  and 
wQ,  will  be  considered.     Hence,  the  discussion  will  refer  to 
the  root  sensitivity  of  the  system.     (Hereafter,  root  sensitivity 
and  pole  sensitivity  will  be  used  interchangeably.) 

Restricting  attention  to  root  sensitivity  does  not  completely 
preclude  the  determination  or  discussion  of  a  best  design.    In 
fact,  there  are  two  arguments  which  partially  justify  (and  might 
even  suggest)  this  simplification.     One  argument  is  that  root 
specifications  are  more  stringent  than  either  the  bandwidth  or 
error  coefficient  specifications.    The  latter  two  are  frequently 
specified  in  terms  of  inequalities,  i.e.  ,  greater  or  smaller  than 
a  certain  minimum  or  maximum  value.     Therefore,   some  tolerance 
can  be  built  into  the  design  by  working  on  the  "safe  side"  of  the 
specifications.    Another  argument  is  based  on  the  root- 
coefficient  relations  of  the  system  characteristic  equation.    If 
the  mobility  of  two  of  the  systems'  closed  loop  poles  are 
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restricted  by  means  of  a  sensitivity  constraint,  then  the  mobility 
of  the  remaining  closed  loop  poles  is  also  restricted.    This 
statement  stems  from  the  fact  that  the  specified  poles  have 
corresponding  roots  which  must  be  factors  of  the  characteristic 
equation.    Now,  in  order  for  these  roots  to  be  perfect  factors  of 
the  characteristic  equation,  the  coefficients  of  the  characteristic 
equation  must  have  remained  reasonably  constant  for  small 
parameter  changes.    This  implies  that  the  remaining  roots,  and 
hence  the  remaining  poles,  also  remained  reasonably  constant. 
Since  both  bandwidth  and  the  error  coefficients  are  partially 
dependent  (and,  in  fact,  totally  dependent  for  low  order  systems 
without  zeroes)  on  the  characteristic  equation,  it  follows  that 
their  sensitivity  functions  bear  an  approximate  relation  to  the 
root  sensitivity  functions. 

Root  Sensitivity  Equations 
Although  the  practical  importance  of  sensitivity  has  long 
been  recognized,  it  has  only  been  in  recent  years  that  effort 
was  expended  in  the  development  of  sensitivity  as  a  useful  tool. 
The  work  of  Kokotovic  and  Siljak,   (2) ,  is  considered  to  be  a 
major  advancement  towards  this  goal.    In  essence,  they  used 
the  definition  of  sensitivity,  equation  (4-1),  in  conjunction  with 
the  generalized  Mitrovic  equations,  equations  (2-6)  and  (2-7), 
to  obtain  algebraic  expressions  for  (what  the  author  considers 
to  be)  the  components  of  root  sensitivity.    Their  results  are  not 
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only  particularly  suited  for  use  in  the  design  of  control  systems 
based  on  a  specified  pole  pair,  but  also  simplify  sensitivity 
computations  because  they  permit  the  exclusive  use  of  real 
variables  in  calculating  the  sensitivity  of  a  function  containing 
complex  variables. 

Since  the  Kokotovic  and  Siljak  sensitivity  equations  will 
be  used  in  subsequent  examples,  a  brief  review  of  their  develop- 
ment will  now  be  presented.    Differentiating  the  generalized 
Mitrovic  equations,  equations  (2-6)  and  (2-7),  with  respect  to 
a  parameter  Q  gives 

OJm      ^"   «j=^  «-'  J    <3      ££  K    , 

+  £$?»4+**irf*o  (4"2) 

+  ^  4s?  *>£  4><  (f)  =  o 

where  f61{  f   )  and  2^-1    '        are  the  derivatives  of  jz£  (  f  )  and 
^k-l(  J^  )  respectively  and  are  related  by  the  general  equation 

4K  (?) = -z[s  &  +*-,$  +■  t^  <C*  wj         (4.4) 

Numerical  values  for  jzf,"  for  various  values  of  f    are  given  in 
Appendix  II.     Manipulation  of  equations  (4-2)  and  (4-3)  and 
application  of  equation  (4-1)  gives 

(4-5) 
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where 


/%?£<**&*  4K<r)     (4- 


Bi  =  £^wf  ^  (f ;  e2=j£  <?«  w$  <?«^  j 


C  -f  & "*  ^.^  ^  -£&■*«««> 


^ 


£Q 


7) 
(4-8) 
(4-9) 


Solving  equations  (4-5)  and  (4-6)  for  SQ  and  S^n  gives 


where 


-    1 

\                               ■ 

if 

D 

(4-10) 

5q     ^~<3 

D 

(4-1 1L 

0   * 

4, 

5; 

a* 

Be, 

(4-12) 

c, 

^2. 

|                            (4-13) 

A, 

A 

E           C* 

(4-14) 

Equations  (4-10)  and  (4-11)  are  the  desired  sensitivity  design 
equations.     (Hereafter,  Sq  and  SqI\  will  be  referred  to  as  the 
"component  sensitivity  functions.  ")    At  this  time,  their  true 
simplicity  is  masked  by  their  dependence  on  the  solution  of  a 
number  of  auxilliary  equations .    Suffice  to  state  that  the  solutions 
of  the  auxilliary  equations  become  exceedingly  simple  when  the 
location  of  a  pole  pair  is  specified  (  f  Q  and  wQ)  and  the 
coefficients  of  the  characteristic  equation  are  known  (aP) .    This 
will  be  clearly  demonstrated  in  subsequent  examples. 
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Intrepretation  of  Sensitivity  Values 

Equations  (4-10)  and  (4-11)  may  be  used  in  two  ways. 
Substitution  of  a  given  set  of  parameter  values  into  the  right 
hand  side  of  the  equations  gives  a  definite  sensitivity  value 
for  a  specific  set  of  parameter  values.    Alternatively,  the 
assignment  of  a  sensitivity  value  to  the  component  sensitivity 
functions  give  the  parameter  relations  which  must  be  satisfied  in 
order  to  obtain  a  system  of  the  specified  sensitivity. 

In  either  case,  the  final  design  will  be  dependent  upon  the 
accepted/assigned  value  to  sensitivity.    It  follows  that  some 
guidelines  for  reasonable  sensitivity  values  must  be  established. 
Since  the  literature  is  lacking  in  this  area,  the  following  guides 
(which  are  consistent  with  the  original  interpretation  of  equation 
(4-1))  are  proposed  by  the  author: 

A.     Signs  (+  or  -)  of  sensitivity  values  indicate  sense 
only.  Therefore,  the  minimum  value  of  sensitivity  is  zero, 
implying  that  a  function  is  totally  insensitive  to  variations  of  a 
parameter. 

B. 

in  Q.i.e.  ,  a  5%  change  in  Q  results  in  a  5%  change  in  P.     This 
is  considered  to  be  the  maximum  permissible  value  of  sensitivity; 
if  possible,   systems  should  have  sensitivity  values  much  less 
than  one. 
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SJJ=  1  indicates  a  linear  variation  in  P  due  to  changes 


c. 


Q 


>1  indicates  amplification  of  the  parameter 
variation.    This  is  generally  considered  undesireable  for  a 
control  system  (but  may  have  application  in  some  detection 
schemes) . 

D.    As  indicated  in  A  above,  sensitivity  values  have  vector 
properties.    Therefore,  the  total  root  sensitivity,  Si^",  is 
appropriately  expressed  as 

(S<RT)2  =  (SF)2+  (Sj^n)2  (4-15) 

Sensitivity  as  a  Specification 

The  concept  of  sensitivity  has  been  developed  to  the  extent 
that  it  can  now  be  used  as  an  additional  system  constraint 
subject  to  analysis  by  the  algebraic  design  method.    The 
remainder  of  this  section  will  be  devoted  to  the  design  and 
discussion  of  specific  systems  which  must  satisfy  sensitivity, 
as  well  as  other,  specifications. 

Example  4 . 1 

Reconsider  example  3.3  in  view  of  a  specification  requiring 
minimum  sensitivity  of  the  specified  poles  with  respect  to 
variations  in  system  gain,  K. 

Considering  the  system  gain,  K,  as  a  parameter,  the  system 
characteristic  equation  is 

(1) 

The  sensitivity  equations  will  now  be  devloped  in  detail  to 

demonstrate  the  solution  of  the  auxilliary  equations.    Applying 

equations  (4-7)  and  (4-8)  gives 

3D 


/?/    =     /  SOO  Oj  —  /OO Oz 

/?2   =  /oo /z az-^t%  o,  —  ?se> i/2  o3 

Bz.  m  /oo  oz  —  2  oco 

after  substituting  wn  =  wQ  =  5(2)1/2  and  values  for  jz£  (0.707) 
and  sz^_^(0.  707)  from  Appendices  I  and  II.     Although  specific 
values  of  P  and  K^  are  unknown,  inspection  of  equation  (1)  in 
conjunction  with  the  results  of  example  3.3  show  that  the 
coefficients  of  the  characteristic  equation  must  have  the 
following  values  for  all  solutions 

QZ=Z£0  Q,     =/££>  /,  <*£  *  1 

which,  when  substituted  into  equation  (2),  gives 

(3) 

Taking  the  partial  of  each  coefficient  in  equation  (1)  with 
respect  to  K  gives 

Applying  equations  (4-9)  and  substituting,  as  above,  f or  wn 
and  jzfk (0.707)  gives 

which,  for  Z  =  5  as  obtained  in  example  3.3,  reduces  to 

c,  *  s- ( / - / o  «r )  Cz  *  $  >"  l£J£T  ~ 0  ^ 

Applying  equation  (4-12)  and  substituting  equation  (3)  gives 

D  »  -ZSo  x  to3  (6) 

Applying  equation  (4-13)  and  substituting  from  equations  (3) 

and  (5)  gives 
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D**1   =  ZSoo  (7) 

Applying  equation  (4-14)  and  substituting  from  equations  (3) 
and  (5)  gives 

£>k    **-'*-&*/& ('-'*#*)  (8) 

Substituting  Q  =  K  =  50,  /    =    fQ  =  0.  707  and  equations (6) , 
(7)  and  (8)  into  equations  (4-10)  and  (4-11)  gives 

^    -M/**r-/J  (9) 

■3**-   Vz  do) 

Equations  (9)  and  (10)  are  the  desired  design  equations. 
Since  only  one  additional  constraint  equation  was  needed,  they 
would  normally  be  squared  and  added  to  obtain  an  expression 
for  the  total  root  sensitivity.    In  this  case,  however,  S^nis 
observed  to  be  a  constant.    Therefore,  the  minimum  pole 
sensitivity  with  respect  to  K  can  be  found  by  setting  equation 
(10)  equal  to  zero.    This  gives  Kj.  =  0.1  which  corresponds  to 
the  solution  of  tachometer  feedback  only. 

For  Kj.  =  0.1,  Swn  =  1/2  and  Sk  =  0.    Therefore,  the  analysis 
predicts  that  a  ten  per  cent  change  in  K  should  produce  a  five 
per  cent  change  in  wR  while  f  remains  constant.    These 
predictions  can  be  mathematically  verified.    The  system 
characteristic  equation  for  the  specified  value  of  Kj.  is 

ForK=50:  wQ=7.07  f    Q=  0.707 

ForK=55:  w    =7.42  Jf  =0.708 
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Therefore  Au>^    ^^9<r%  Aj_  =   0m/+% 


Parameter  Selection 
The  preceeding  example  raises  several  questions  which 
have  yet  to  be  discussed.    The  first  deals  with  the  selection  of 
a  parameter  on  which  to  base  the  sensitivity  analysis.    One 
choice  might  be  the  parameter  which  is  considered  to  have,  the 
widest  variations  during  expected  conditions  of  operation,  e.g.  , 
the  gain  of  an  electronic  amplifier.    Another  choice  might  be 
similar  to  the  first  but  restricted  to  plant  parameters.    This 
assumes  the  designer  is  able  to  select  and/or  control  the 
sensitivity  of  the  required  compensating  elements.    There  is  no 
guarantee  that  either  is  correct.    It  may  be  possible  that  the 
system  is  naturally  insensitive  to  the  parameter  likely  to  have 
the  widest  variations  in  value  while  it  could  be  highly  sensitive 
to  a  parameter  that  is  expected  to  have  only  minute  variations 
in  value.    It  is  also  possible  that,  by  minimizing  the  sensitivity 
due  to  one  parameter,  the  effect  of  another  is  >  amplified*    The 
answer,  then,  requires  investigation  of  all  possible  parameters. 
It  is  believed  that  a  thorough  investigation  of  a  particular  type 
of  system  could  lead  to  the  establishment  of  general  guides  for 
the  selection  of  the  appropriate  parameter(s)  in  other  similar 
systems.    Since  such  an  investigation  is  not  consistent  with  the 
purpose  of  this  paper,  the  author  will  be  content  to  arbitrarily 
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select  any  parameter  which  will  facilitate  demonstration  of  the 
algebraic  method. 

It  should  also  be  noted  that  there  are  no  restrictions  as  to 
the  number  of  parameters  which  may  be  similtaneously  considered. 
If  the  sensitivity  of  pole  P  with  respect  to  parameters  Q^ ,  Q2 
and  Q3  is  desired,  the  component  sensitivity    functions  for 
each  parameter  may  be  individually  calculated .    The  total 
component  sensitivity  function  due  to  all  three  parameters  is 
the  algebraic  sum  of  the  individual  components,  e.g.  , 

.    Since  the  sense 
of  each  component  must  be  considered,  there  is  a  possibility 
that  individual  components  may  cancel.    Another  fortunate 
circumstance  is  that  the  computations  for  a  situation  considering 
the  root  sensitivity  with  respect  to  N  parameters  do  not  increase 
N-fold.    This  is  because  equations  (4-7),   (4-8)  and  (4-12),  i.e., 
Aj,  A2,  B]_,  B2  and  D,  remain  the  same  for  all  parameters. 

Minimum  Sensitivity 
Another  point  of  interest  in  example  4.1  is  the  ease  with 
which  the  condition  of  minimum  sensitivity  was  obtained.    This 
is  not  a  unique  situation  because  it  was  possible  to  obtain 
results  by  inspection  in  other  systems  studied  also.     However, 
it  seems  quite  likely  that  for  many  systems  the  parameter 
relations  for  minimum  sensitivity  must  be  obtained  by  setting 
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the  derivative  of  the  sensitivity  function  equal  to  zero.     (Note 
the  problem  if  the  sensitivity  with  respect  to  more  than  one 
parameter  is  being  considered.)    In  such  cases,  care  must  be 
taken  to  insure  that  the  function  to  be  differentiated  is  in  a 
form  which  will  indicate  minima  near  zero.     This  is  necessary 
because  minimum  sensitivity  was  defined  to  be  zero.    The  above 
stipulation  can  be  conveniently  accomplished  by  considering 
either  the  absolute  value  or  the  square  of  the  component 
sensitivity  function,  i.e.,    /S^^j         or   (Sq n )      vice    S<?  ** 
Use  of  the  total  sensitivity  function,  equation  (4-15),  is  also 
appropriate. 

Discussions  thus  far  have  repeatedly  used  the  term  "minimum 
sensitivity."    This  does  not  necessarily  mean  that  minimum 
sensitivity  must  be  attained  by,  or  even  an  objective  of,  every 
design.    The  specified  sensitivity  value,  whether  it  be  zero  or 
two,  should  be  consistent  with  the  other  specifications.  For 
example,  if  the  pole  locations  are  only  required  to  be  accurate 
within  five  per  cent  for  satisfactory  system  response  and  the 
parameter  of  interest  is  also  expected  to  have  five  per  cent 
variations,  a  sensitivity  value  of  one  is  sufficient.    Of  course, 
if  analysis  shows  that  a  lower  value  is  attainable,  the  design 
should  be  based  on  the  lower  value  provided  it  can  be  achieved 
without  additional  cost  or  degradation  of  some  other  aspect  of 
system  performance. 
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The  nature  of  the  sensitivity  function  in  the  vicinity  of  the 
design  value  is  also  of  concern.    If  the  design  value  is 
associated  with  a  point  on  the  sensitivity  function  curve  where 
the  curve  is  very  steep  or  irregular,  the  system  is  likely  to  be 
more  sensitive  than  an  alternative  design  based  on  a  (reasonably) 
higher  value  which  corresponds  to  a  smooth  portion  of  the 
sensitivity  function  curve. 

Example  4.2 

A  portion  of  a  control  system  has  the  forward  transfer 
function,  G(s)  =  K/s(s  +  2)  (s  +4).    Dominant  closed  loop  poles 
with  jT0=  1/2  and  w0  =  5  are  required  for  the  desired  response. 
Overshoot  is  of  prime  importance;  therefore ,     C  Q  should  be 
maintained  within  1%  accuracy  for    [v  anticipated  gain  variations 
of  5%.    Velocity  and  acceleration  signals  are  required  elsewhere 
in  the  system  and  are  therefore  available  for  use  as  feedback 
signals. 

The  problem  as  stated  gives  3  specifications  (pole  pair 
and  SjJ  and  3  parameters.     (Note  that  an  additional  specification 
concerning  Swn  would  normally  require  the  insertion  of  an 
additional  parameter  whereas  a  specification  regarding  the  total 
pole  sensitivity  would  not.    Also  note  that  the  specifications 
permit  wQ  to  vary,  if  necessary,  in  order  to  maintain  f  Q 
approximately  constant.) 
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System  CE:      s3  *•  2*  (k  +  <K<x}  ts(ftKkT)t  K  =  © 
Mitrovic  equations:  K  *    zs  *"  zs  K  Ka  (1) 

Applying  equations  (4-7)  thru  (4-9): 


(2) 


(3) 


At   -  31S-  Sao*.  /? a.  -  SOqz. -  So, 

C^i-ZSKKv  C*,  =S(£Ka-Kr)       J 

Unlike  example  4.1,  the  coefficients  of  the  characteristic 
equation  do  not  have  fixed  numerical  values.    Therefore,  the 
solution  of  the  determinants  for  D  and  D,  ,  equations  (4-12)  and 
(4-14),  will  be  hampered  by  the  presence  of  all  3  parameters. 
At  this  point  it  will  be  convenient  to  solve  equations  (1)  and  (2) 
for  K    and  K.  .    Then  equation  (3)  can  be  expressed  in  terms  of  K 
only. 

«n  I2y  Kr=     ^#£  (4) 

Therefore 

/),  =  IZS-2K  Az    =-   K  ■*-  IZS 

C,    ^  25/k  Cz  *  -I/o/k 

7 
Using  equation  (5)  to  solve  for  D  and  0^: 

Using  equation  (6)  to  solve  equation  (4-13: 


^<  K*  —  IZSK  ?  I2SZ 


(7) 


Inspection  of  equation  (7)  shows  that  s£  =  0  for  K  =  86.5.     Since 
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the  only  requirement  is  that  S,*  $:  0.2,  this  gain  setting  should 
be  considered  as  a  possible  design  value.      For  K  =  86.5, 
Kj.  =  0 .  397  and  KQ  =  0 .  0284 .    For  these  values ,  the  third  closed 
loop  pole  is  located  at   d"  =  -3.4.    Therefore,  the  complex  poles 
are  not  completely  dominant.    In  addition,  figure  4-1  shows  the 
sensitivity  function  is  very  steep  at  the  point  of  minimum 
sensitivity.    Again  using  figure  4-1,  the  system  will  be  designed 
on  the  basis  of  K  =  600.    The  other  design  values  are:  Kj-  =  0.238 
and  Ka  =  0.0384.    The  complex  poles  are  assured  of  being 
dominant  because  the  third  pole  is  now  located  at  <T~=  -24 .    It 


w„ 


is  stated  without  proof  that  S^n  also  approaches  zero  as  K 

f 

becomes  very  large.    It  is  also  less  sensitive  than  S/ ,  having 

a  value  <    0.05  at  the  design  gain  setting. 
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*5£/js/7-/  vry  (S$)      vs      G,*/n    (k) 
/-or.       £x^^^/.s         4--£ 
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5.      Conclusions 

The  purpose  of  this  paper  was  to  develop  a  design  technique 
based  on  the  premise  that  design  can  be  accomplished  by  solving 
a  set  of  algebraic  simultaneous  equations  derived  from  the 
system  design  specifications.    In  addition  to  demonstrating  the 
feasibility  of  this  concept,  the  examples  in  the  preceeding 
sections  also  showed  that  the  method  developed  has  certain 
desirable  characteristics.    Specifically,  the  algebraic  method 
is  analytically  simple,  capable  of  handling  any  number  of 
parameters,  rapid,  specific  and  accurate. 

The  method  is  in  no  way  restricted  to  the  exclusive  use  of 
the  design  specifications  considered  in  this  paper.    The 
specifications  considered  were  simply  selected  to  enhance  the 
presentation.    In  general,  a  system  can  be  designed  to  satisfy 
any  set  of  independent  specifications  providing  the  specifications 
can  be  expressed  as  algebraic  equations.    It  is  obvious,  of 
course,  that  the  solution  of  the  resulting  equations  becomes 
more  difficult  as  the  number  of  specifications  increase.    For  this 
reason,  it  is  suggested  that  analysis  be  normally  limited  to 
specifications  which  are  considered  to  be  prime,  i.e.  ,  those 
which  are  exactly  specified  as  opposed  to  others  stated  in  terms 
of  inequalities  or  approximations.    The  results  of  this  analysis 
can  then  be  used  to  determine  if  the  secondary  specifications  are 
also  satisfied. 
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The  inclusion  of  sensitivity  equations  into  the  set  of  system 
constraint  equations  was  considered  especially  significant.    Here 
again,  the  use  of  sensitivity  equations  is  in  no  way  dependent 
upon  the  assumptions  made  in  Section  4  concerning  the  relation 
of  root  sensitivity  to  system  sensitivity,  i.e.  ,  use  of  root 
sensitivity  as  a  performance  index.    This  was  logically 
demonstrated  in  example  4 . 2  where  a  root  sensitivity  component 
was  of  particular  concern  while  the  system  sensitivity  had  no 
restrictions  whatsoever.    Although  only  root  sensitivity  was 
discussed,  there  is  no  apparent  reason  why  similar  restrictions 
cannot  be  imposed  on  other  design  specifications. 

As  previously  indicated  in  Section  4,  there  is  a  definite  need 
for  the  study  of  specification/parameter  sensitivity  relationships 
before  sensitivity  can  become  a  useful  design  tool.    At  the 
present,  the  sensitivity  effects  of  all  possible  parameters  must 
be  evaluated  with  respect  to  each  specification  and  to  each  other 
before  a  meaningful  design  can  be  achieved.    Analysis  of  even 
a  simple  system  on  this  basis  would  require  a  prohibitive  amount 
of  computation  because  of  the  large  possible  number  of  parameter/ 
specification  combinations  available. 

"Seek  Simplicity  and  Distrust  It" 
Alfred  North  Whitehead 
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